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We adopt a continuous model to estimate the Grothendieck constants. An analytical formula to 
compute the lower bounds of Grothendieck constants has been explicitly derived for arbitrary orders, 
which improves previous bounds. Moreover, our lower bound of the Grothendieck constant of order 
three gives a refined bound of the threshold value for the nonlocality of the two-qubit Werner states. 

PACS numbers: 03.65.Ud, 03.67.-a 


I. INTRODUCTION 

Quantum mechanics exhibits the nonlocality of the 
nature in essence. The impossibility of reproducing all 
correlations observed in composite quantum systems us¬ 
ing models a la Einstein-Podolsky-Rosen (EPR) [l] was 
proven in 1964 by Bell Q. A quantum state is said to 
admit a local hidden variable (LHV) model if all the mea¬ 
surement outcomes can be modeled as a classical random 
distribution over a probability space. Consider a bipar¬ 
tite state p in Hu ® H% with subsystems A and B. If Alice 
performs a measurement A on the subsystem A with an 
outcome oq and, at space-like separation, Bob performs a 
measurement B on the subsystem B with an outcome bj , 
then an LHV model supposes that the joint probability 
of getting di and bj satisfies 

Pr (ai,bj\A, B, p) = I Pr(a,i\A, \)Pr(bj\B, \) du p (\), 
Jn 

where dui p { A) is some distribution over a space H of hid¬ 
den variable A. A quantum state is called local if it admits 
an LHV model, and nonlocal otherwise. 

In his seminal work, Bell showed that all quantum 
states admitting LHV models satisfy the so-called Bell 
inequalities [2]. That is, a state admits no LHV mod¬ 
els if it violates some Bell inequalities. It is known that 
every pure entangled bipartite or multipartite state vio¬ 
lates a generalized Bell inequality SI- Namely, for pure 
states the entanglement and the non-locality coincide. 
However for mixed states, the situation is more compli¬ 
cated. There are no general methods to judge whether a 
mixed state admits an LHV model or not, i.e. to find all 
Bell inequalities is computationally hard [a, @]. Even for 
the most concerning two-qubit Werner states, the precise 
threshold value of nonlocality is still unknown. 

As the “fundamental theorem in the metric theory of 
tensor products”, the Grothendieck’s theorem Q had a 
major impact on Banach space theory. The constants 
related to the Grothendieck’s theorem are nowadays 
called Grothendieck constants @]- It turns out that the 


Grothendieck constants are related to the Bell inequal¬ 
ities, observed by Tsirelson Of particular interest, 
Acin-Gisin-Toner 0 demonstrated that the threshold 
value for the non-locality of the Werner states for projec¬ 
tive measurements is given explicitly by the Grothendieck 
constant of order three, K( 3), see Section [n] for the def¬ 
inition. This reduces the problem of the nonlocality of 
Werner states to the computation of the exact value of 
K{ 3). However, it is formidably difficult to compute the 
Grothendieck constants except for the case of order two. 
Generally, what one can do is to estimate the lower and 
upper bounds of the Grothendieck constants. 

In this paper, by generalizing the Bell OTierator with 
465 measurement settings on each side in [1 lj] to a contin¬ 
uous model with infinitely many measurement settings, 
we present an analytical formula in estimating the lower 
bounds of the Grothendieck constants. This formula is 
valid for Grothendieck constants of arbitrary order and 
improves many previously obtained bounds. From our 
lower bound of K( 3), we derive a bound of the threshold 
value for the nonlocality of the two-qubit Werner states, 
which gives the best knowledge about such nonlocality 
up to date. 


II. LOWER BOUND OF GROTHENDIECK 
CONSTANT FOR ARBITRARY ORDER d 


Let A4 m (R) be the set ofraxm real matrices and 5 d_1 
the unit sphere in m, d £ N. Given M £ Af m (R), we 
define 


C\M ) = sup 


^ ^ - 1 1 : rj aibj , 

*>l=i 


(1) 


where the supremum is taken over all possible assignment 
a *> bj G {1,-1}, 1 < i,j < m. Replacing ai,bj by el- 
dimensional unit vectors, we define 


Q(M) = sup 


.bieS 1 *- 1 


Mij • bj 

*,l=i 
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where the supremum is taken over all d-dimensional unit 
vectors a; and bj and a, • bj denotes their scalar product. 
The Grothendieck constant of order d is defined by 


K(d) = sup sup 

m> 1 M£Mm 
M jtO 


Q(M) 
(R) C(M) 


( 2 ) 


The value of C(M) depends only on the choice of the 
matrix M. Besides M, Q{M) also depends on the di¬ 
mension d of Euclidean space where we choose unit 
vectors ai and bj. It is a great challenge to evaluate 
the Grothendieck constant K(d) for general d. Till now 
the only exactly known result of K(d) is for d = 2, 
K( 2) = y/2 0. For d > 3, there are some lower 
bounds of K(d): For instance, Briet-Buhrman-Toner 0 
obtained a lower bound of K(d) for general d, 


K (d) > ^ 
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In this section, we propose a continuous model to 
compute the lower bounds of the Grothendieck constant 
for arbitrary d. Our results improve some known lower 
bounds. In particular, the lower bound of K{ 3) we ob¬ 
tained is the best up to date, which improves the result 
on the non-locality of the Werner states. 

For any n G N, let m = n + n(n — l)/2. We choose the 
following special M € M m ( R) such that 0 


C{M) = sup 


n 

£ a-ibj + Y^[ a ij( b i 

i,j=l i<j 


bj) + o,j)] , 


where the supremum is taken over a*, aij G {1,-1} and 
bj,ftij € {1,-1} ({di,aij} stands for {aj}™ 1 in the def¬ 
inition ©, similarly {6,;, #3^} for {&j}JT 1 ). Hence by 
choosing atj = sign(&i — bj) and /3jj = sign(ai — aj), 
one has 


C(M) = sup 


£ a i b s+J2(\ b i 

i,j=l i<j 


bj | + | di 



Using this discrete version of optimization problem, 
Vertesi [0 runs a numerical simulation to obtain the 
lower bounds of K(d), 3 < d < 5. Since the complexity of 
the computation grows exponentially, the computation of 
the lower bounds becomes impossible for large d. Instead 
of this, we propose a continuous model to reformulate this 
problem. Let P(S d ~ 1 ) denote the space of probability 
measures on Then we have 


Theorem II. 1. For d > 1, 

K(d) > sup^p^-i) (l/sd-rXdMx) 


+ / |x 

Js d ~ 1 xS d ~ 1 


y| dn(x)dji( y) . 


(5) 


Proof. From the discrete version 0, the assertion 0 
holds for rational convex combination of delta measures, 
i.e. /x = with Aj G Q+, = 1 

and ai G £ d_1 . Note that any probability measure 
G P(S d _1 ) can be approximated by convex combina¬ 
tion of delta measures in the weak topology (precisely, 
weak* topology in the terminology of functional analy¬ 
sis). This follows from Krein-Milman theorem and the 
fact that delta measures are extreme points for the set of 
probability measures (13]. Since |x|, |x — y| are continu¬ 
ous functions, the theorem follows from the definition of 
weak convergence of measures. □ 


Hence, to derive an effective lower bound of K (d) it suf¬ 
fices to choose some good measures for this optimization 
problem. It seems that the problem becomes more com¬ 
plicated since the finite dimensional optimization prob¬ 
lem has been transformed to an infinite dimensional prob¬ 
lem on the space of measures. However, this is in fact an 
advantage which allows one to choose nice absolutely con¬ 
tinuous measures on a sphere to get explicit lower bounds 
of the Grothendieck constants. 

Let (0i, 02, • • • , 4>d- 1) be the spherical coordinates of 
S"* -1 C such that 

( Xi = sin0i sin0 2 • • • sin(0 d _i), 
x 2 = sin0i sin0 2 • • • cos(0d_i), 


Correspondingly, for given M, to get a better bound 
of Q(M ) one needs to suitably choose the vectors a, and 
bj. By setting vectors ay (/3y resp.) parallel to bi — bj 
(ai — aj resp.) and ai = bi, 1 < i < n, one obtains 


Q(M)= sup^ ( |^ ai | 2 +2^|ai-aj| ] . (3) 


i<j 


Therefore 


K (d) > sup 
ai eS d ~ 


£ a > 


+ T2£l ai_a jl ’ ( 4 ) 


i^j 


Xd-i = sin0i cos0 2 , 
l Xd = cos 01, 

where (pi G [0,7r] for 1 < i < d — 2 

and <pd -1 G [0, 27t). We denote by dji = 

sin d ~ 2 0i sin d-3 02 • • • sin tpd -2 d<p\d<p 2 ■ ■ ■ depd -i the 

spherical (volume) measure of S'® -1 . For simplicity, we 
also denote by x = (x',Xd), x' G R d_1 , the Cartesian 
coordinates of x G S d_1 , by Pd(x) = Xd the projection 
to the d-th coordinate, and by 0i(x) the first spherical 
coordinate of x. For any 0 < a < b < §,we denote 

!l‘ = {x G 5 ,d_1 : a < 0i(x) < b}, 


where the supremum is taken over (ai, • • • ,a n ) G 
(S ,d_1 ) n for any n > 1. 


which is contained in the upper hemisphere S d 1 D {xd > 
0}. To obtain the lower bounds of K(d ), we will choose 










3 


the uniform probability measure on 17^, i.e. fi a ,b = 
voin b v °l|n b ■ The variables a, & are introduced to refine 
(i.e. maximize in some sense) the lower bound of K(d), 
since we have no priori knowledge of the optimal mea¬ 
sure which attains the maximum on the right hand side 
of (0. 

Then second term of the right hand side of 0 corre¬ 
sponds to 


ln b „ xn b 


|x - y\dfx a: b(x)d^a,b(y) 


vol(fl^) 2 J n b xn b 


:-y|d/r(x)d^(y), 


which involves a 2 (d— l)-multiple integration. This is the 
obstruction for the numerical computation of the lower 
bounds for large d. Nevertheless, the computation can be 
considerably simplified by the symmetry of the sphere. 

Lemma 1. For any x, x £ S ,d_1 satisfying </>i(x) = 
</>i(x), we have 


Since SO(d) acts isometrically on S d 3 , the equation © 
follows. □ 


By the Lemma, for given x £ the integral over y 
only depends on ^i(x). Without loss of generality, we 
may choose x = ( 0 , 0 , ••• , sin </>i, cos i.e. its spher¬ 
ical coordinate is {<j> i, 0, - - - ,0). Then for any y whose 
spherical coordinate reads (V’i, ip2, ■ • ■ , ifd- i), one has 


|x — y| = y/2 — 2 sin fa sin ipi cos ^2 — 2 cos <f> i cos ifi. 

(7) 

Therefore we obtain 


In* xn b 


|x - y\dn(x)dn{y) 


= vol(S ' 2 3 )vol(S d 3 ) f f f /(^i,V’i,V’ 2 )#idV’i #2 
J a J a J 0 



l x — y|d/r(y) 



|x-y|d/r(y). 


( 6 ) 

for d > 3, where 


Proof. For given x and x with P^(x) = Pd(x), there is an 
isometry of S"* -1 , A £ SO(d ), such that H(x) = x and 
A(fik) = 17^. This can be seen as follows: Without loss 
of generality, one may write x = (x', c) and x = (x', c) 
where c = P^x). Since x',x' £ y/l — c 2 S d ~ 2 , there is a 
B £ SO(d — 1) such that Bx' = x'. Hence one can take 



/((/>!, ^ 1 ,^ 2 ) = |x-y| sin d 2 f> 1 sin 6 * 2 V’i sin d 3 if 2 , ( 8 ) 

with |x — y| given by ©• Combining above results, we 
have the following theorem: 

Theorem II.2. The Grothendieck constant K(d), d> 3, 
satisfies 


K{d) > 


(/ a & sin d 2 (fidtfi^j 


sin d 1 (j) 1 


d- 1 


vol(5' 


d-3\ r b r b r n 


vo\(S d ~ 


m i r 

/ ( 0 i, ^ 1 , ipi)d(f>i #1 #2 


(9) 


where /(</>i, ^> 2 ) is given by 0 and ©. 


By reducing a 2 (d— l)-multiple integration to a triple 
integration for any d > 3, we have obtained a lower 
bound of the Grothendieck constants which can be eas¬ 
ily calculated via numerical methods. By varying a, b 
in the domain {(a, b) : 0 < a < b < £}, one may 
get a refined lower bound of K(d) by maximizing the 
right hand of © for these a, b. For instance, by taking 
a = 0, b = 1.04819755, one gets K{ 3) > 1.41758. Taking 
a = 0.742832, b = 0.749115, one gets K( 5) > 1.46112. 
Some numerical results are listed in the following table: 



K{d) > 

d 

Vertesi 

Briet etc. 

our result 

3 

1.41724 

1.33333 

1.41758 

4 

1.44521 

1.38791 

1.44566 

5 

1.46007 

1.42222 

1.46112 

6 


1.44574 

1.47017 

7 


1.46286 

1.47583 

8 


1.47586 

1.47972 

9 


1.48608 

1.48254 


The first column collects the results of Vertesi 0 for 
the case d = 3,4,5. The second column contains the 
lower bounds of K(d) proved by Briet-Buhrman-Toner 
0 - From this table, one immediately fi gur es out that 
our results improve the results of Vertesi [ll| and Briet- 
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Buhrman-Toner 0 for d < 8. 


III. NON-LOCALITY OF TWO-QUBIT 
WERNER STATES 


In 1989, Werner explicitly constructed LHV models for 
some entangled mixed bipartite states 0] . The two-qubit 
Werner state is given by 


P^=P\i/j )(ip |+(l-p)I/4, 


( 10 ) 


where in computational basis, \ip~) = (|01) — |10))/-\/2 
and 0 < p < 1. is separable if and only if p < 1/3 
[0. It admits an LHV model for all measurements for 
p < 5/12 0, and admits an LHV model for projective 
measurements for p < 0.6595 [13- 
Let A t and Bi, i = 1,2, be dichotomic observ¬ 

ables with respect to the two qubits, Ai = a i ■ a and 
Bi = b, ■ a, with a = (oi, U 2 , 03) the Pauli matri¬ 
ces. a i = (af^af^af 5 ), b, = {b { ~ i L \bf\bf ) ) are 3- 
dimensional real unit vectors. For any Bell operator, 


B{M) = MijAiftBj, 
*,j=i 


where M G A4 m (K) as in ©, the mean value is given by 

m 

Tr(B(M)pW)=pJ2Mi J aj-bj. 
i,j =1 

Therefore the maximal violation of the corresponding 
Bell inequality is given by p K (3) and p™ admits LHV 
models for projective measurements if and only if p < 
i/m) 0. Hence the nonlocality problem of the two- 
qubit Werner states is reduced to estimate the value of 
A(3). 

However, the precise value of K{ 3) is still unknown. 
There are various attempts to derive the upper and 
lower bounds of the Grothendieck constants. For 
instance, Krivine 0 ] showed that K{ 3) < 1.5163. 
The Clauser-Horne-Shimony-Holt (CHSH) inequality im¬ 
plies that I\( 3) > y/2, TJ. Vertesi 0 constructed 
Bell inequalities involving 465 settings on each qubit 
m > 465) to show that K{ 3) > 1.417241, 
i.e., p^p admits no LHV models for p > 0.705596. From 
Theorem 111.21 we have shown that K{ 3) > 1.41758. 
Therefore p^ admit no LHV models for p > 0.705428. 
This provides the best known bound for the nonlocality 
of two-qubit Werner states, (see Fig. 1: Acin-Gisin- 
Toner show that p/) admits LHV models for p < 0.66. 
Vertesi shows that p™ admits no LHV models for p > 


0.705596. We show that p™ admits no LHV models for 
p > 0.705428.) . . .. 


Acin-Gisin-Toner 


Our work 

Wemer 



Vertesi 

Barrett 


CHSH 

Separable 


->• 


1 5 

3 12 

i °' 66 f 

0.705428 / 

\ 0.705596 


FIG. 1: Nonlocal properties of two-qubit Werner states. 


IV. CONCLUSION AND DISCUSSIONS 

We have presented an analytical formula to estimate 
the lower bounds of the Grothendieck constants for arbi¬ 
trary order. It has been shown that our bounds improve 
the previously obtained bounds for d = 3,..., 8. It is also 
straightforward to calculate the lower bounds for higher 
order d. However, the lower bounds of K{d) for large d 
derived by this approach cannot exceed 1.5, which has 
been proven in 11|, section IV. This certainly doesn’t 
beat the best lower bound of 1.677 for K( 00) in 0 . 

Nevertheless, our lower bound of K( 3) gives a bound 
of the best threshold value for the nonlocality of the two- 
qubit Werner states up to date. In fact, our new lower 
bounds of the Grothendieck constant of high orders can 
be also used to improve the knowledge about the non¬ 
locality for higher dimensional quantum states such that 
the related mean values of Bell operators are determined 
by the Grothendieck constants 0 . 

So far there are no effective ways to justify whether 
a quantum entangled state admits local hidden vari¬ 
able models or not. For two-qubit Werner states, the 
Grothendieck constant K (3) plays the essential role that 
p™ admits LHV models if and only if p < K( 3) -1 . It 
would be also interesting to find such quantities for gen¬ 
eral quantum states. 
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